non-ergodic mesoscopic systems 
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Suppose there is a mesoscopic system connected to single channel leads. If the system is non- 
chaotic or non-ergodic then the thermodynamic and transport properties do not depend on impurity 
averaged density of states. We show that the partial density of states as well as density of states of 
a given system can be determined exactly from the asymptotic wave-function (or scattering matrix) 
at the resonances. The asymptotic wave-function can be determined experimentally without any 
knowledge about the quantum mechanical potential (including electron-electron interaction) or wave 
function in the interior of the system. Some counter intuitive relations derived here can allow this. 
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For bulk samples, ensamble averaging makes it un- 
necessary to know the exact impurity configuration and 
the exact Hamiltonian of the sample. This is because ex- 
periments can only observe ensamble averaged physical 
properties. Ensamble averaging works because of the well 
known ergodicity hypothesis and the fact that sample to 
sample fluctuation is not very large [![. For mesoscopic 
systems even when dealing with an ensamble, sample to 
sample fluctuations are often so large that we cannot talk 
of any averaged physical quantity [2| . If disorder averag- 
ing cannot be taken then it becomes necessary to know 
the exact impurity configuration of a system, a seemingly 
impossible task for an experimentalist. Also sometimes 
for mesoscopic systems ergodicity itself does not hold for 
a single sample that is further illustrated below. In this 
work we show some exact relations. This relation is not 
only counter intuitive but can also provide a way to ex- 
perimentalists, of bypassing the knowledge of the exact 
Hamiltonian of the system to arrive at thermodynamic 
and transport properties of mesoscopic system. 

Consider an arbitrary potential V(x,y) in the shaded 
region £1 of Fig. 1. The z degree of freedom is usually 
frozen due to strong confinement in the z direction. This 
potential can define a mesoscopic system that could be 
a quantum dot or a quantum ring or anything else ||. 
Typically such a system is coupled to leads or measuring 
probes. The leads connect the system to reservoirs that 
are at fixed chemical potentials. The reservoirs inject (or 
absorb) electrons to (or from) the system through the 
leads. We consider single channel leads as most experi- 
ments are done with single channel leads Q • When leads 
are multichannel, then we do not have enough control 
over the properties of the system to make it of practical 
use 0. By taking specific forms of V(x, y) we can make 
the system to leads coupling to be weak or strong. For 
a fully chaotic (or ergodic) system, electrons will access 
all the states in the system specially when the coupling 
to the leads is very weak and the electron spends suffi- 
cient time in the system. However, for a non-chaotic (or 
non-ergodic) system, all the states will not be accessed. 
Only part of the states will be accessed that will de- 
pend on the position and the details of the leads and the 



system, that will constitute the partial density of states 
(PDOS). So the contribution of these electrons to ther- 
modynamic observables like quantum capacitance of the 
system, or heat capacity of the system, as well as the lin- 
ear response non-equilibrium effects will be determined 
by this PDOS. This PDOS cannot be determined from 
the Hamiltonian of the isolated system as the PDOS de- 
pends on initial conditions (that is through which lead 
the electrons enter the system and through which lead 
they leave as well as the characteristics of the leads and 
the system). However, the scattering matrix depends on 
these factors. So the scattering matrix has more infor- 
mation than the Hamiltonian and scattering matrix for- 
mulation is very important for mesoscopic systems Q. 

The many body Hamiltonian for a N particle system 
is a function of 3N coordinates. It has been proved that 
the motion of one of these N particles is governed by an 
effective potential that is just a one body potential and 
a function of 3 coordinates [j| • So in the preceding para- 
graph, when we refer to the confining potential V(x, y) of 
the system, we are referring to this one body effective po- 
tential that includes electron-electron interaction exactly. 
The determination of this one body potential is however 
very difficult and has never been exact up to date 
For bulk systems there are fairly good although approx- 
imate methods for obtaining the effective potential, but 
the same need not be applicable to finite systems. So for 
mesoscopic systems if we can bypass the determination 
of the internal details and the exact Hamiltonian of the 
system, by using the S-matrix then the S-matrix includes 
the effect of electron-electron interaction exactly. 

The approach proposed in this paper is due to some re- 
cent experiments 0, H| that are motivated by the possibil- 
ity of obtaining important information from the scatter- 
ing phase shift. So the present work is an effort to identify 
information that can be obtained from such experiments. 
A series of experiments [9|, [lfj also tell us that resonances 
in such a system as that schematically shown in Fig. 1 
are Fano resonances. Recently it has been shown for 
some particular potential (namely a delta function po- 
tential) in a quantum wire, that at the Fano resonance 
the density of states (DOS) and partial density of states 
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(PDOS) can be determined exactly from semi-classical 
formulas involving scattering phase shift, although Fano 
resonance is a purely quantum phenomenon [ll| . In this 
work we show that these results are true for any general 
potential V(x,y). 

We shall use the following identity to arrive at our con- 
clusions. It can be rigorously derived for one dimension 
(ID) that is naturally also true for quasi one dimension 
(Q1D) [l2j that involves tracing over sub-band index. 
In simpler form, [l3| for any such system schematically 
shown in Fig. 1. 
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Here A denotes scattering matrix and E m denotes the 
kinetic energy in the mth channel that is further defined 
below. 

The Schrodinger Eq. describing the system is 
h 2 d 2 d 2 

' " (2) 
For \x\ < a, V(x,y) is negative enough to create at least 
one bound state (otherwise there will be no Fano reso- 
nance). For \x\ > a, V(x,y) — V c (y), where V c (y) is 
the confinement potential in the leads. The most general 
solution to Eq. 2 in different regions of Fig. 1 is 
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Here 6 n (y) is the solution to the following Eq. that holds 
for |x| > a. 



~^ + V c (y)]e n (y) = e n 9 n (y) 



Therefore, 
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If e n > E then one can see from Eq. 6 that k n is imag- 
inary (the mode is evanescent). The electron is incident 
along the mth channel which implies e m < E. For all 
n < to, e„ < E and such channels are propagating. For 
n > (to + l),e„ > E, implying that these channels are 
evanescent. For n < to, r nm and t nm gives reflection and 
transmission amplitudes, respectively. They also consti- 
tute the elements of the scattering matrix A. For n > to, 
t nm and r nm gives transition amplitude from the prop- 
agating mode to the nth evanescent mode in the right 
lead and left lead, respectively. However, they are not 
scattering matrix elements. 



We can now define two functions ^ ev \x,y) and 

^ od \x,y) such that ^ ev \x,y) = ip( ev \-x,y) and 
tl) {od ^{x,y) = -^ od \~x,y). 

^ ev \x, y) = ^ =1 (S nm e~^ x - Sffie**"*)^ (7) 

V "ti 

^ od \x, y) = ^ =1 (S nm e' lk ^ - S^e ik - x ) 6 -^r (8) 
Then one can see that as well as is given by 
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This works because any function can be written as a sum 
of an even function and an odd function. And any square 
matrix can be written as a sum of a symmetric matrix 
and an antisymmetric matrix. 

Due to the same principle, the wave-function in the 
scattering region f2 can be written as a sum of an even 
function and an odd function. We denote them as 
4>n V \x,y) and 4>n d \x,y). 

4>^\x,y) = ^ =l c mX t v) {x,y) (12) 

^° d \x,y) = V™ =l c mX t d) {x,y) (13) 

where x^ ev ^ and ^( od ) satisfy the Schrodinger Eq. 2 with 
V(x, y) = in the region fi. Therefore an( l X 

will satisfy the following Eqs. 



xtHa,y) = 0m(y) for \y\<b 
X ^°\x,y)=0fory = F(x) or G(x) 

x<£ v) (x,y) = xt ) (-x,y) 



xt d) {x,y) = -xT ) i-x,y) 



(14) 
(15) 
(16) 
(17) 



F(x) and G(x) define the two curves at the upper and 
lower boundaries of the region fl. Here 'eo' stands for 'ev 
or od\ We also define the following matrix elements 
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Now we require (j)n°^ and 



dx 



to be continuous at 



x = a for all \y\ < b. Thus we get 

y^oo is —ik m a c(eo) ^k rn a\ m 
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= iC m 8 m (y) 
(19) 
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Multiplying Eq. 19 and Eq. 20 with \6 n (y) and inte- 
grating from y = — b to y = b and then combining them 
we get a single matrix equation. 
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or 



S<>£) = e" lfc - a [l + 2«(F( e °) - il)- 1 ] mn e- lk " a (22) 

It is known in scattering theory that the bound states 
of the potential V(x, y) can also be obtained from Eqs. 
7 and 8 by omitting the terms 5 nm e tknX . Without this 
term Eqs. 7 and 8 are solutions to the Sc. Eq. in 2 with 
correct boundary condition wherein there is no incident 
wave. Identical analysis that lead to Eq. 21 in this case 
gives 
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Here mt is the threshold value of m for which bound 
states exist. And n m = ik m . For m < mt states will 
be scattering states. Supposing only the first channel is 
propagating then mt = 2. This will be further illustrated 
soon. Solutions to Eq. 23 or solutions to the following 
Eq. will give bound states. 
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Here 'cc' means closed channel. Let us partition f ( eo ) 
into propagating and evanescent (or closed) channels. 



Therefore, 
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where 
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If there is only one propagating channel then m = n = 
1 . Also for one propagating channel, p = 1 . Hence G^ 
becomes a number. Therefore from Eq. 28 
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All the infinite S mn appear in Su through F cc , which 
carry information of the entire Hilbert space accessed by 
the incident wave. For non-ergodic systems there can be 
many more states not accessed by the incident wave, ie., 
the transition amplitude between these states and the 
incident wave being 0. These states are not populated. 
So they neither contribute to DOS and PDOS nor to 
the scattering matrix. Trivial examples are cases like 
when the incident wave has a certain symmetry that is 
incompatible with the symmetry of these states. More 
generally, there can be localized and scarred states that 
are known to occupy a subset of the Hilbert space and 
do not connect to the rest. 

Threshold energy E for 1st closed channel is given by 
£2) > 0. Below this energy the 2nd channel can 
have bound states. Such bound states will occur at en- 
ergies given by the solution to Eq. 24. At these energies 
the 1st channel will be propagating as its threshold is 
given by jjjt{E — ei) > and Su is given by Eq. 30. But 
at these energies G eo will diverge as it includes matrix 
elements of [F c e c ° — il]^ 1 as can be seen from Eqs. 31 and 
24. That in turn implies that at a Fano resonance (as 
can be seen from Eq. 32) 
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From this one can show that 
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So from Eq. 22 and Eq. 27 (for m and n being propa- 
gating channels) 
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This is consistent with the fact that at Fano resonance 
tn=0 (which is the definition of Fano resonance) as can 
be seen from Eqs. 11 and 33. Therefore tu — tu=0 
at Fano resonance. Also from Eqs 10 and 33 at Fano 
resonance 



m - r* n = -i[sin{25 {e) ) + sin{25 {o) 
Therefore, at Fano resonance 
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Therefore at Fano resonance, from Eqs. 1 and 35 
1 
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Here Arg(A mn ) — Arctan I ^^ mn The complicated inte- 
gral involving the local potential on LHS is the PDOS 
for an electron incident in channel n and scattered to 
channel m. So at the Fano resonance the PDOS can 
be determined exactly from the scattering phase shift 
d \- Ar 9(^'^)] and |A mn | 2 . Both these quantities were mea- 
sured in 0, H [l(J. Summing over m and n we get at 
Fano resonance 



Am mn J n mn SV(r) 



HC = 



d[^-logDet[A]] 
dE 

(37) 

which is Friedel sum rule. It is expected to hold good 
only for bulk samples that are in semiclassical limit. For 
finite systems that are in quantum regime there is always 
a correction term arising from j^—(A mn — A^ n ). We 
do not know of any example where this correction term 
becomes exactly 0. In Q1D, although Fano resonance is a 
purely quantum interference phenomenon, the correction 
terms are exactly making Friedel sum rule exact. The 
correction terms are extremely non- universal (resonances 
are generally characterized by line shape i.e., |A mn | 2 and 
scattering phase shifts Arg(A mn )) and also depend on 
sample specific parameters (e.g., E m depends on V c (y) 
as well as material parameters like bottom of conduction 
band, effective mass fi, etc.) 



In the experiments of Refs. @, d, G3 ~jEArg(t) turns 
out to be negative at the Fano resonance. This means 
that PDOS at a Fano resonance is negative. That in 
turn means that the quantum capacitance at a Fano reso- 
nance can be negative. Which means the effective poten- 
tial (V = Q/C) due to some negatively charged electron 
in a quantum dot at a Fano resonance can be positive. 
So negatively charged electrons in the leads will be at- 
tracted or in other words there will be electron-electron 
attraction at a Fano resonance of a quantum dot. Such 
an attraction was observed in numerical simulations [l4| , 
although a proper explanation could not be given. 

As a simple experiment one can repeat the experiments 
of Ref. @, H, EH with distorted quantum dots for which 
the spacing between the peaks in |A mn | 2 are not expected 
to be uniform. This spacing will depend on DOS which 
for distorted dots do not peak at regular intervals. Thus 
from the spacings one can get the LHS of Eq. 36. Also 
from similar data as that in Refs. @, d, Gi| one can get 
RHS of Eq. 36 as explained before. 

In conclusion, Friedel sum rule and similar semiclassi- 
cal formulas (Eqs. 36 and 37) become exact at resonances 
for any general potential in single channel Q1D that can 
support a resonance. This is very counter intuitive as 
Fano resonance is a purely quantum phenomenon. We 
do not know of any other situation where semiclassical 
formula can become exact as mesoscopic world is always 
quantum. Semiclassical formula can at most be a good 
approximation. Thereby, the experimental data of Refs. 
0i S E3] do carry important information. Using such 
data experimentalists can bypass the knowledge of mi- 
croscopic Hamiltonian to know the DOS and PDOS and 
hence thermodynamic and transport properties. 
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